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ABSTRACT 

We compare the peculiar velocities measured in the SBF Survey of Galaxy Distances with the pre- 
dictions from the density fields of the IRAS 1.2 Jy flux-limited redshift survey and the Optical Redshift 
Survey (ORS) to derive simultaneous constraints on the Hubble constant Hq and the density parameter 
P = n°-^/b, where b is the linear bias. We find /3/ = 0.42^°;^° and /3o = 0.26±0.08 for the /i?^5'and ORS 
comparisons, respectively, and ifp = 74 ± 4 km s~^ Mpc (with an additional 9% uncertainty due to the 
Cepheids themselves). The match between predicted and observed peculiar velocities is good for these 
values of Hq and /?, and although there is covariance between the two parameters, our results clearly 
point toward low-density cosmologies. Thus, the unresolved discrepancy between the "velocity-velocity" 
and "density-density" measurements of /3 continues. 

Subject headings: galaxies: distances and redshifts — cosmology: observations - distance scale - 
large-scale structure of universe 



1. INTRODUCTION 

Gravitational instability theory posits that the present- 
day peculiar velocity Vp of a galaxy will equal the time- 
integral of the gravitational acceleration due to nearby 
mass concentrations, i.e., Vp = {g)t. For galaxies out- 
side virialized structures, linear perturbation theory can 
be used to relate Vp to the present-day local mass density: 
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(Peebles 1980), where ft is the mass density of the universe 
in units of the critical density, Hq is the Hubble constant, 
and Smir) is the mass density fluctuation field. A fur- 
ther common simplification is the linear biasing model, 
Sg = bSm, where 6g is the fluctuation field of the observed 
galaxy distribution, and b is the linear bias. With these 
assumptions, the observed peculiar velocity will be pro- 
portional to the quantity (3 = D,°-^/b. 

Redshift surveys of complete samples of galaxies can 
be used to determine the galaxy density in redshift space 
Sg{z), which is then smoothed and used to predict the 
peculiar velocity field. In this case, the distances are mea- 
sured in units of the Hubble velocity and Hq cancels out 
of Eq. (1). Comparison with observed pecuhar velocities 
from distance surveys tied to the Hubble flow then yields 
a value for f3 (see Strauss & Willick [1995] for a compre- 
hensive review of the methods) . Recent applications using 
TuUy-Fisher distances and the density field of the IRAS 
1.2 Jy redshift survey (Fisher et al. 1995) have found best- 
fit values I3i = 0.4-0.6 (Schlegel 1995; Davis et al. 1996; 



1998), where the subscript denotes the IRAS survey. Riess 
et al. (1997) have done the comparison with Type la su- 
pernovae (SNIa) distances and find /3/ — 0.40 ± 0.15. On 
the other hand, the Potent method (Dekel et al. 1990), 
based on the derivative of Eq. (1) but incorporating nonlin- 
ear terms, gives a /?/ about twice this value, most recently 
Pi = 0.89 ± 0.12 (Sigad et al. 1998). 

The surface brightness fluctuation (SBF) method 
(Tonry & Schneider 1988) of measuring early-type galaxy 
distances is new to the field of cosmic flows. A full review 
of the method is given by Blakeslee et al. (1999). The 
/-band SBF Survey (Tonry et al. 1997, hereafter SBF-I) 
gives distances to about 300 galaxies out to ~ 4000 km s~^. 
Tonry et al. (1999, hereafter SBF-II) used these data to 
construct a parametric flow model of the surveyed volume, 
but as the analysis made no use of galaxy distribution in- 
formation, it could not constrain fl. It did, however, allow 
for a direct tie between the Cepheid-calibrated SBF dis- 
tances and the unperturbed Hubble flow, and thus a value 
for Hq. Previous Hq estimates with SBF (e.g., Tonry 
1991; SBF-I) relied on far-field ties to methods such as 
TuUy-Fisher, Dn-(J, and SNIa, and thus incurred addi- 
tional systematic uncertainty. Unfortunately, this "direct 
tie" of SBF to the Hubble flow was still not unambiguous. 
For instance, Ferrarese et al. (1999, hereafter F99) used 
a subset of the same SBF data but a different parametric 
flow model to derive an Hq 10% lower than in SBF-II. 

The present work uses SBF Survey distances in an ini- 
tial comparison to the peculiar velocity predictions from 
redshift surveys and derives simultaneous constraints on 
Hq and f3. The comparison also provides tests of linear 
biasing, the SBF distances, and claimed bulk flows. 
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2. ANALYSIS 

We compute the redshift-distance (cz-d) relation in the 
direction of each sample galaxy using the method of Davis 
& Nusser (1994), which performs a spherical harmonic so- 
lution of the redshift-space version of the Poisson equation. 
Computations are done for both the IRAS 1.2 Jy (Fisher 
et al.l995) and "Optical Redshift Survey" (Santiago et al. 
1995) catalogs as a function of their respective density pa- 
rameters Pi and f3o- The smoothing scale of the cz-d 
predictions is typically '^ 500 km s^^. 

The predictions arc then compared to the SBF obser- 
vations using the x^ minimization approach adopted by 
Riess et al. (1997), who assumed a fixed "redshift error" 
of (T„ = 200 km s~^. This term includes uncertainty in 
the linear predictions, unmodeled nonlinear motions, and 
true velocity measurement error; we also report results 
with cr„ = 150 km s~^. We use the same subset of galax- 
ies as in SBF-II, namely those with high-quality data, 
(V—I)o > 0.9, and not extremely anomalous in their ve- 
locities (e.g., Cen-45); we also omit Local Group members. 
The full sample then comprises 280 galaxies. Unlike the 
case for the SNIa distances, we lack a secure external tie 
to the Hubble flow, so we successively rescale the SBF 
distances to different values of Hq and repeat the mini- 
mization. 

The Nusser-Davis method cannot reproduce the mul- 
tivalued cz-d zones of clusters, but it can produce ma- 
jor stall regions in the cz-d relation. Since such regions 
will have high velocity dispersions, this deficiency in the 
method may be overcome with allowance for extra velocity 
error. A significant fraction of SBF galaxies reside in the 
Virgo and Fornax clusters. We adopt the following three 
approaches in dealing with these. "Trial 1" uses individual 
galaxy velocities but allows extra variance in quadrature 
for the clusters according to: aci{r) = a^j \J\ -f- (r/ro)^, 
where o-q = 700 (400) km s"^ and r^ = 2 (1) Mpc are 
adopted for Virgo (Fornax). These spatial profiles are 
meant to mimic the observed projected profiles and param- 
eters from Fadda et al. (1996) and Girardi et al. (1998); 
at large-radius ctci follows the infall velocity profile (e.g., 
SBF-II). "Trial 2" uses a fixed velocity error but removes 
the virial dispersions by assigning galaxies their group- 
averaged velocities. The groups are defined in SBF-I: 29 
galaxies are grouped into Virgo, 27 into Fornax, 7 into 
Ursa Major, and all other groups have 2-6 members, with 
34% of the sample being ungrouped. "Trial 3" is similar 
to trial 2, but 72 sample galaxies within 10 Mpc of Virgo 
(about 50% larger than the zero-velocity radius found in 
SBF-II) and 30 within 5 Mpc of Fornax are removed, re- 
ducing the sample by 36% to 178 galaxies. 

For each of the trials, we calculate x^ for the compar- 
isons with both the IRAS and ORS gravity fields, using 
(T„ of both 150 and 200 km s~^, for a range in Hq in steps 
of 0.1 km s~^ Mpc and a range in /3 in steps of 0.1. We 
then interpolate in f3 by cubic splines to find the Hq-/] 
combination minimizing x^- 

3. RESULTS 

Figure shows the 68%, 90%, and 99% joint confi- 
dence contours on Hq and /3 from the x^ analyses for 
trial 2 (which gives the median best-fit Hq) using ay = 200 
km s~^. Significant covariance exists between the two pa- 



rameters. Table 1 shows the x^ minimization results for 
all the different comparison runs. In order to give a more 
realistic impression of the true uncertainties, the tabulated 
errors are for (5x^ — 2.3, the 68% joint confidence range 
on 2 parameters, except for the last column which varies 
Hq within uncertainty limits to explicitly take account of 
the covariance in /3 (see below). 
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Fig. 1.— The 68%, 90%, and 99% joint probability contours on 

t) and f3 = QP-^ /b are plotted for the "Trial 2" comparison (see 
for description) of the SBF peculiar velocities against the IRAS 
(solid curves at right) and ORS (dotted curves at left) density field 
predictions. Crosses mark the best-fit positions. 

Table 1 
Results from Peculiar Velocity Comparisons 



Run 


CTi, 


y2 


/9 


± 


Ho 


± 


feo = 


=74 ± 


IRAS-1 . 


200 


0.99 


0.39 


0.06 


73.2 


1.2 


0.42 


+0.10 
-0.07 


IRAS-1 . 


150 


1.21 


0.40 


0.07 


73.7 


1.1 


0.41 


+0.09 
-0.06 


IRAS-2 . 


. 200 


0.89 


0.43 


0.06 


73.9 


1.0 


0.43 


+0.10 
-0.06 


IRAS-2 . 


150 


1.18 


0.45 


0.06 


74.3 


1.0 


0.44 


+0.08 
-0.05 


IRAS-3 . 


. 200 


1.05 


0.40 


0.08 


74.1 


1.2 


0.40 


+0.10 
-0.05 


IRAS-3 . 


150 


1.33 


0.41 


0.07 


74.4 


1.2 


0.41 


+0.07 
-0.05 


ORS-1.. 


. 200 


1.05 


0.23 


0.05 


72.1 


1.2 


0.25 


+0.07 
-0.05 


ORS-1.. 


. 150 


1.28 


0.24 


0.05 


72.4 


1.0 


0.27 


+0.06 
-0.05 


ORS-2.. 


. 200 


0.97 


0.25 


0.05 


72.8 


1.3 


0.29 


+0.06 
-0.07 


ORS-2.. 


. 150 


1.29 


0.27 


0.05 


73.1 


1.2 


0.30 


+0.05 
-0.05 


ORS-3.. 


. 200 


1.12 


0.22 


0.05 


74.1 


1.4 


0.21 


+0.05 
-0.04 


ORS-3.. 


. 150 


1.43 


0.24 


0.06 


74.6 


1.5 


0.22 


+0.05 
-0.04 



Columns list: the comparison run (designated as survey— trial) ; 
assumed small-scale velocity error (T„ (km s~^); reduced x^ 
(degrees of freedom = number of galaxies minus 2); best-fit 
/3 and Hg (km s"! Mpc); and best-fit /3 for //q = 74 ± 1.4. 

On average, the IRAS comparisons prefer /3/ « 0.4 and 
Hq « 74, while the ORS, which more densely samples the 
clusters, prefers f3o ~ 0.25 and Hq « 73. Of course, there 
can be only one value of Hq for the sample galaxies. The 
table shows that when the cluster galaxies are removed in 
trial 3, the best-fit Hq increases to 74 for the ORS, while 
it changes little for IRAS. We adopt if o = 74 ± 1.4 as the 
likely value from this velocity analysis, where we use the 
median Hq error for the IRAS trials combined in quadra- 
ture with the variance in Hq among the trials. The error 
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increases to ±4 km s^^ Mpc when the 5% statistical uncer- 
tainty in the tie to the Cepheid distance scale is added in 
quadrature (SBF-II). Finally, the estimated uncertainty in 
the Cepheid scale itself is ~9% (F99; Mould et al. 1999). 
These uncertainties in Hq due to the distance zero point 
have no effect on the uncertainty in /3. 
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Fig. 2. — The reduced x^ is plotted as a function of /3 for ii'o = 74 
km s~^ Mpc and cr„ = 200 km s~^. Circles connected by solid spline 
curves show results from "Trial 2," which includes all 280 galax- 
ies, and diamonds connected by dashed curves show results from 
"Trial 3," which removes 102 galaxies near the Virgo and Fornax 
clusters. Solid/open symbols are used for the /iJAS/ORS compar- 



Figure || illustrates how the reduced x^ for Hq — 74 and 
CTi, — 200 varies with /3 for trials 2 and 3 (spanning the 
range in best-fit /3 and x^) of the IRAS and ORS com- 
parisons. Overall, the last columns of Table 1 indicate 




+600 km/s 



Resid 



Pi 



0.42 



-1-0.10 
-0.06 



and Po = 0.26 ± 0.08, where the errors 



include the uncertainty for a given trial (including the 2% 
uncertainty in Hq from the velocity tie) and the variation 
among the trials. The results are independent of whether 
one takes a median or average of the trials, or adopts either 
ay — 150 or 200 km s~^. However, /3/ and /?o would in- 
crease for example by ~ 30% to 0.56 and 0.33, respectively, 
if Hq were known to be 5% larger for a fixed distance zero 
point, or decrease by 15-20% if Hq were 5% smaller, but 
again, the x^ analysis indicates Hq is constrained to 2% 
for a fixed distance zero point. 

Finally, Figure ^ compares the /3/ = 0.4 predicted and 
Hq = 74 observed peculiar velocities in the Local Group 
frame for all the galaxies, using their group-averaged ve- 
locities to limit noise. Given that the two sets of peculiar 
velocities were derived independently (apart from our ad- 
justment of Hq and /9/), the agreement is quite good. The 
predictions resemble a heavily smoothed version of the ob- 
servations, and xt confirms this, although the residuals 
may hint at a slightly misaligned dipole. Large negative 
residuals near {I, b) ^ (283°, 4-74°) indicate Virgo backside 
infall in excess of what the spherical harmonic solution can 
produce. These issues will be addressed in more detail by 
a forthcoming paper (Willick et al. 2000, in preparation). 



Fig. 3. — Predicted peculiar velocities from the IRAS survey 
(top) with [3j = 0.4, observed peculiar velocities from the SBF sur- 
vey (middle) with Hq = 74 km s~^ Mpc, and the residuals (observed 
minus predicted, bottom) are shown in Galactic coordinates and the 
Local Group frame. Light/dark circles show negative/positive pe- 
culiar velocities, coded in size as shown; smaller symbols are plotted 
on top of larger ones for clarity. Unlike the predictions, the obser- 
vations have no smoothing applied and contain both distance and 
velocity errors, causing the greater range in symbol size. 



4. DISCUSSION 

We have found the most consistent results for Hq ~ 74 
km s~^Mpc, intermediate between the favored values of 
77 and 71 reported by SBF-II and the Hq Key Project's 
analysis of SBF (F99), respectively. The differences in 
these "SBF Hq" values result entirely from the flow mod- 
els^. SBF-II explored a number of parametric flow models 
ranging from a pure Hubble flow plus dipole to models 
with two massive attractors and a residual quadrupole, 
as well as some with a local void component. The pure 
Hubble flow model gave the same Hq as obtained by F99. 
Adding Virgo and Great Attractors gave Hq — 73.5, and 
the additional Local Group-centered quadrupole gave a 
further 6% increase in Hq. Each component significantly 
improved the model likelihood, but it was suggested that 
the quadrupole arose from inadequate modeling of the flat- 
tened Virgo potential. If so, then Hq may be better esti- 
mated using a sharply cutoff, Virgo-centered quadrupole. 
Figure 23 of SBF-II shows that such a model would indeed 
yield iJo « 74. 

^We have changed the Key Project value by 1.6% to be appropriate for our distance zero point, since our concern here is with the tie to the 
Hubble flow. 
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We also note that most of the SBF-II models had 
a significant excess Local Group peculiar motion of '^ 
190 km s~^, unlike the IRAS density field predictions (e.g., 
Willick etal 1997). However, when the excess was modeled 
as a "push" from a large nearby void, Hq dropped from 
78 to 73 because of the underdensity introduced. This 
model achieved the best likelihood of any considered, but 
the treatment of the void was deemed too ad hoc to qual- 
ify as a standard component of the flow model. Thus, it 
may be that the SBF-II result for Hq suffered from the 
arbitrariness inherent in parametric modeling. However, 
the good match we find between the predicted and ob- 
served velocity fields supports the claim by SBF-II that, 
after accounting for the attractor infalls, any bulk flow of 
the volume cz < 3000 km s^^ is < 200 km s^^, as the 
structure within this volume accounts for most of our mo- 
tion in the cosmic microwave background rest frame (see 
Nusser & Davis 1994). 

The best-fit values of /3/ == 0.421°^° and Po = 0.26 ± 
0.08 are similar to the (3i = 0.40 ± 0.15 and f3o = 0.3 ± 0.1 
results found by Riess et al. (1997) using 24 SNIa dis- 
tances out to cz ~ 9000 km s"^, nearly 3 times our survey 
limit. It is interesting that despite the close agreement 
on /3, SBF and SNIa still disagree by nearly 10% on Hq 
(e.g., Gibson et al. 1999). This implies that the discrep- 
ancy is mainly due to the respective distance calibrations 
against the Cepheids (otherwise the Hq offset would cause 
major disagreement on /3). We also note that the ratio 
f^i/Po ~ 1.6 agrees well with the estimates by Baker et al. 
(1998). 



Our results are consistent with all other recent compar- 
isons of the gravity and velocity fields (so-called "velocity- 
velocity" comparisons), regardless of the distance sur- 
vey used. Thus, the near factor-of-two discrepancy with 
the Potent analysis (a "density-density" comparison) of 
Sigad et al. (1998) persists. Since the latter analysis is 
done at larger smoothing scales, one obvious explanation 
is non-trivial, scale-dependent biasing, and some simula- 
tions may give the needed factor-of-two change in bias over 
the relevant range of scales (Kauffmann et al. 1997; but 
see Jenkins et al. 1998). However, at this point, the incon- 
sistency in the results of the two types of analysis remains 
unexplained. 

We are currently in the process of analyzing the SBF 
data set using methods that deal directly with multival- 
ued redshift zones and incorporate corrections to linear 
gravitational instability theory to take advantage of SBF's 
ability to probe the small-scale nonlinear regime (Willick 
et al. 2000, in preparation). We also plan in the near fu- 
ture to use the data for a "density-density" determination 
of (3. In addition, we are working towards a direct, unam- 
biguous tie of SBF to the far-field Hubble flow, which will 
significantly reduce the systematic uncertainty in /3. 
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